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ABSTRACT 



A diversity of 2:1 resonance configurations can be expected in extrasolar planetary 
systems, and their geometry can provide information about the origin of the resonances. 
Assembly during planet formation by the differential migration of planets due to planet- 
disk interaction is one scenario for the origin of mean-motion resonances in extrasolar 
planetary systems. The stable 2:1 resonance configurations that can be reached by 
differential migration of planets with constant masses and initially coplanar and nearly 
circular orbits are (1) anti-symmetric configurations with the mean- motion resonance 
variables 6i = \i — 2A2 + wi and 62 = \i — 2A2 + ^2 (where \j and Wj are the mean 
longitudes and the longitudes of periapse) librating about 0° and 180°, respectively (as in 
the lo-Europa pair), (2) symmetric configurations with both 61 and 62 librating about 0° 
(as in the GJ 876 system), and (3) asymmetric configurations with 9i and O2 librating 
about angles far from either 0° or 180°. There are, however, stable 2:1 resonance 
configurations with symmetric {61 k, O2 ~ 0°), asymmetric, and anti-symmetric (^i ~ 
180° and 62 ~ 0°) librations that cannot be reached by differential migration of planets 
with constant masses and initially coplanar and nearly circular orbits. If real systems 
with these configurations are ever found, their origin would require (1) a change in 
the planetary mass ratio mi/m2 during migration, (2) a migration scenario involving 
inclination resonances, or (3) multiple-planet scattering in crowded planetary systems. 
We find that the asymmetric configurations with large 62 and the Oi ~ 180° and 62 ~ 0° 
configurations have intersecting orbits and that the 61 ^ 62 ^ 0° configurations with 
ei > 0.714 have prograde periapse precessions. 



Mean-motion resonances may be ubiquitous in extrasolar planetary systems. There are as 
many as three resonant pairs of planets in the 12 multiple planet systems discovered to date. It is 
well established that the two planets about the star GJ 876 discovered by Marcy et al. (2001) are in 
2:1 orbital resonance, with both lowest order, eccentricity-type mean-motion resonance variables. 
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(where Ai,2 are the mean longitudes of the inner and outer planet, respectively, and wi^2 are the 
longitudes of periapse), librating about 0° (Laughlin & Chambers 2001; Rivera & Lissauer 2001; Lee 
& Pcale 2002). The simultaneous librations of 9i and 02, and hence the secular apsidal resonance 
variable 

= Wi - W2 = 01 - 02, (3) 

about 0° mean that the periapses are nearly aligned and that conjunctions of the planets occur 
when both planets are near periapse. The other two pairs of planets suspected to be in mean- 
motion resonances are the pair about HD 82943 in 2:1 resonance (Mayor et al. 2004) and the inner 
two planets about 55 Cnc in 3:1 resonance (Marcy et al. 2002). 

The symmetric geometry of the 2:1 resonances in the GJ 876 system was not expected, because 
the familiar 2:1 resonances between the Jovian satellites To and Europa are anti-symmetric, with 
01 librating about 0° but 02 and ^3 librating about 180° (a configuration that would persist in the 
absence of the Laplace resonance with Ganymede). In the lo-Europa case, the periapses arc nearly 
anti-aligned, and conjunctions occur when lo is near periapse and Europa is near apoapse. Lee 
&; Peale (2002) have shown that the differences in the resonance configurations are mainly due to 
the magnitudes of the orbital eccentricities involved (see also Beauge & Michtchenko 2003) . When 
the eccentricities are small, the lo-Europa configuration is expected from the resonant perturbation 
theory to the lowest order in the eccentricities. However, the eccentricities of the GJ 876 system are 
sufficiently large that there are large contributions from higher order terms, and stable simultaneous 
librations of 0i and 02 r6C[u.ir6 0]^ ^ O2 ^ 0° for a system with eccentricities and masses like those 
in GJ 876. 

A scenario for the origin of mean-motion resonances in extrasolar planetary systems is that they 
were assembled during planet formation by the differential migration of planets due to gravitational 
interaction with the circumstellar disk from which the planets formed. A single giant planet (with 
a planet-to-star mass ratio > a few xl0~^) can open an annular gap in the circumstellar gas disk 
about the planet's orbit. If a disk forms two giant planets that arc not separated too far (with 
the ratio of the orbital scmimajor axes ai/a2 ~ 1/2), the planets can also clear the disk material 
between them rather quickly (Bryden et al. 2000; Kley 2000; Kley, Peitz, k Bryden 2004). Disk 
material outside the outer planet exerts torques on the planet that are not opposed by disk material 
on the inside, and the outer planet migrates toward the star. Any disk material left on the inside of 
the inner planet exerts torques on the inner planet that push it away from the star. The timescale 
on which the planets migrate is the disk viscous timescale, whose inverse is (Ward 1997) 



^ = -aiH/afn 
2a2 2 ^ ' ' 



(4) 
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where a = da/dt, the kinematic viscosity v is expressed using the Shakura-Sunyaev a prescription 
[u = aH^Q), H is the scale height of the disk, and J7 = 27r/P and P ^ 2Tra^/^ / {Gmo)^/^ are, 
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respectively, the mean motion and period of an orbit of semimajor axis a. Although the depletion of 
the inner disk means that the inner planet may not move out very far, the condition of approaching 
orbits for capture into mean-motion resonances is established. 

Lee (fe Pcalc (2002) have shown that the observed, symmetric, 2:1 resonance configuration 
in the GJ 876 system can be easily established by the differential migration due to planet-disk 
interaction (see also Snellgrove, Papaloizou, & Nelson 2001). They have also found that the observed 
eccentricities of the GJ 876 system require either significant eccentricity damping from planet-disk 
interaction or resonance capture occurring just before nebula dispersal, because continued migration 
of the planets while locked in the 2:1 resonances can lead to rapid growth of the eccentricities if 
there is no eccentricity damping. As we shall see in §4, there are other types of stable 2:1 resonance 
configurations with large eccentricities for a system with the inner-to-outer planet mass ratio, 
mi/ 1712, of the GJ 876 system, but these configurations are either unstable for planets as massive 
as those in GJ 876 or not reachable by differential migration of planets with constant masses and 
coplanar orbits. 

As Lee & Peale (2002) pointed out, a practical way to investigate what stable resonance con- 
figurations are possible from continued migration of the planets after capture into 2:1 resonances is 
through numerical migration calculations without eccentricity damping. Either eccentricity damp- 
ing or the termination of migration due to nebula dispersal would lead to a system being left 
somewhere along a sequence. They noted that asymmetric resonance configurations with 6i and 
02 librating about angles other than 0° and 180° are possible for systems with mi/m2 different 
from those in GJ 876. Lee & Peale (2003a) have presented some results on the variety of stable 
2:1 resonance configurations that can be reached by differential migration (see also Ferraz-Mello, 
Beauge, & Michtchenko 2003). Asymmetric librations were previously only known to exist for ex- 
terior resonances in the restricted three-body problem (e.g.. Message 1958; Beauge 1994; Malhotra 
1999). Other recent works on 2:1 resonance configurations for two-planet systems include those by 
Beauge, Ferraz-Mello, k Michtchenko (2003), Hadjidemetriou & Psychoyos (2003), Ji et al. (2003), 
and Thommes Sz Lissauer (2003). With the exception of Thommes Sz Lissauer (2003), who exam- 
ined the possibility of inclination resonances for non-coplanar orbits, all of the works mentioned in 
this paragraph have focused on systems with two planets on coplanar orbits. 

In this paper we present the results of a series of direct numerical orbit integrations of planar 
two-planet systems designed to address the following questions. What are the possible planar 
2:1 resonance configurations? Can they all be obtained by differential migration of planets? In 
§2 we describe the numerical methods. In §3 we present the stable 2:1 resonance configurations 
that can be reached by differential migration of planets with constant masses. The analysis in 
§3 is an extension of those by Lee &; Peale (2003a) and Ferraz-Mello et al. (2003). We consider 
configurations with 0.1 < mi/m2 < 10 and the orbital eccentricity of the inner planet, ei, up to 
0.8 and provide a detailed analysis of properties such as the boundaries in mi/m2 for different 
types of evolution, the region in parameter space for configurations with intersecting orbits, and 
the region where the retrograde periapse precessions expected for orbits in 2:1 resonances are 
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reversed to prograde. By using a combination of calculations in which mi / m2 is changed slowly and 
differential migration calculations, we show in §4.1 that there are also symmetric and asymmetric 2:1 
resonance configurations that cannot be reached by differential migration of planets with constant 
masses. In their studies of systems with masses similar to those in HD 82943, Ji et al. (2003) and 
Hadjidemetriou &; Psychoyos (2003) have discovered anti-symmetric 2:1 resonance configurations 
with Oi 180° and ^2 ~ 0°. We show in §4.2 that a sequence of this configuration with anti-aligned, 
intersecting orbits exists for all mi/m2 examined. In §5 wc summarize our results and describe 
several mechanisms by which the configurations found in §4 can be reached. 

2. NUMERICAL METHODS 

Wc consider a system consisting of a central star of mass mo, an inner planet of mass mi, 
and an outer planet of mass m2, with the planets on coplanar orbits. We refer to the planet with 
the smaller orbital semimajor axis as the inner planet, but it should be noted that some of the 
configurations studied below have intersecting orbits and the outer planet can be closer to the star 
than the inner planet some of the time. In addition to the mutual gravitational interactions of the 
star and planets, the direct numerical orbit integrations presented below include one of the following 
effects: forced orbital migration, a change in mi/m2, or a change in mi + m2. We change mi/m2 
according to dln{mi/m2)/dt = constant and keep m^ + m.2 constant for the calculations with a 
change in mi/ 1712, and vice versa for the calculations with a change in mi + m2. For the calculations 
with forced orbital migration (and constant masses), we force the outer planet to migrate inward 
with a migration rate of the form 02/02 oc P2^- The migration rate in equation (4) is of this form if 
a and H/a are independent of a. This form of the migration rate also has the convenient property 
that the evolution of systems with different initial outer semimajor axis 02,0 but the same initial 
01/02 is independent of 02,0 if we express the semimajor axes in units of 02,0 and time in units 
of the initial outer orbital period ^2,0- Except for the migration calculations discussed in the last 
paragraph of §3.2, there is no eccentricity damping. 

The direct numerical orbit integrations were performed using modified versions of the sym- 
plectic integrator SyMBA (Duncan, Levison, & Lee 1998). SyMBA is based on a variant of the 
Wisdom-Holman (1991) method and employs a multiple time step technique to handle close en- 
counters. The latter feature is essential for the integrations that involve intersecting orbits and/or 
become unstable. Lee &: Peale (2002) have described in detail how SyMBA can be modified to 
include forced migration and eccentricity damping. If there is no eccentricity damping, a single 
step of the modified algorithm starts with changing 02 according to the forced migration term for 
half a time step, then evolving the system for a full time step using the original SyMBA algorithm, 
and then changing 02 according to the forced migration term for another half a time step. For 
the calculations with a change in mi/m2 or mi + m2, we modify SyMBA to include the change in 
mi/m2 or mi -|- m2 in the same manner as the forced migration term. 

Unlike the calculations in Lee Sz Peale (2002), which used astrocentric orbital elements, we have 
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input and output in Jacobi orbital elements, and wc apply the forced migration term to the Jacobi 
02 in the differential migration calculations. Lcc & Peale (2003b) have shown that Jacobi elements 
should be used in the analysis of hierarchical planetary systems where 01/0,2 is small and that the 
use of astrocentric elements can introduce significant high-frequency variations in orbital elements 
that should be nearly constant on orbital timescales. We have found through experiments that 
Jacobi elements are useful even for 2:1 resonant systems with 01/02 ~ 2~^/^. In their differential 
migration calculations for the G J 876 system with (mi + 777.2) /'t^o = 0.0095, Lee & Peale (2002) did 
not sec the libration of 02 about 180°, which is expected for small eccentricities. However, when we 
repeated the calculations in Jacobi coordinates, we were able to see the libration of 62 about 180° 
at small eccentricities with the reduced fluctuations in the orbital elements (see also §3.2 below). 

3. RESONANCE CONFIGURATIONS FROM DIFFERENTIAL MIGRATION 

OF PLANETS WITH CONSTANT MASSES 

3.1. Results for (mi + 1712) /mo = 10"^ and d2/a2 = -10-^/^2 

Wc begin with a series of differential migration calculations with 1711/1712 ranging from 0.1 to 
10. The total planetary mass [mi + 1712) /rriQ = 10^'^. The planets arc initially on circular orbits, 
with the ratio of the semimajor axes 01/02 = 1/2 (far from the 2:1 mean-motion commensurability) , 
and the outer planet is forced to migrate inward at the rate 02/02 = — 10'~^/P2- The small total 
planetary mass [compared to (mi -|- 1712) /tuq ^ 10~^ of the GJ 876 system] and slow migration 
rate (compared to eq. [4] with the nominal parameter values) are chosen to reduce the libration 
amplitudes. The slow migration rate is also chosen to reduce the offsets in the centers of libration 
of the resonance variables due to the forced migration. The effects of a migration rate similar to 
that in equation (4) and a larger total planetary mass are discussed in §3.2. 

Figure 1 shows the evolution of the semimajor axes Oj, eccentricities ej, and resonance variables 
9i and 62 for the calculation with mi/m2 = 0.3 (which is close to mi/m2 = 0.32 of the GJ 
876 system). As the outer planet migrates inward, the 2:1 mean- motion commensurability is 
encountered, and the system is first captured into the 2:1 resonance configuration with 0\ librating 
about 0° and 62 about 180° (as in the lo-Europa pair). Continued migration forces ei to larger values 
and 62 from increasing to decreasing until 62 ~ when ei ~ 0.1 (see Fig. 56 for a better view of this 
behavior). Then the system passes smoothly over to the configuration with both 9i and 62 librating 
about 0° (as in the GJ 876 system), and both ei and 62 increase smoothly. The system remains 
stable for e\ up to 0.8, when the calculation was terminated, although it may become unstable 
at larger e\. In Figure 2a we show the anti-symmetric resonance configuration with 6\ 0° and 
62 ~ 180°. The eccentricities of the ellipses representing the orbits are exaggerated so that the 
positions of the periapses (marked by the dashes) are more visible. The planets (represented by 
the small dots) are shown at conjunction. The periapses are nearly anti-aligned, and conjunctions 
occur when the inner planet is near periapse and the outer planet is near apoapse. In Figure 26 
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Fig. 1. — Evolution of the semimajor axes, ai and 02, eccentricities, ei and 62, and 2:1 mean- 
motion resonance variables, 61 = Xi — 2X2 and ^2 = Ai — 2X2 +'^2, for a differential migration 
calculation with mi/m2 = 0.3 and (mi + m2)/mo = 10~^. The outer planet is forced to migrate 
inward with 02/02 = — 10~^/P2- The semimajor axes and time are in units of the initial orbital 
semimajor axis, 02,0, and period, P2,o of the outer planet, respectively. The sequence of resonance 
configurations after resonance capture is (^1,^2) ~ (0°, 180°) — (0°,0°). 



we show the symmetric resonance configuration with 0i ~ ^2 ~ 0°. The eccentricities are for the 
configuration near t/P2,o = 6 x 10^ in Figure 1. The periapses are nearly aligned, and conjunctions 
occur when both planets are near periapse. The evolution shown in Figure 1 is characteristic for 
mi/m2 ^ 0.95. 

Figure 3 shows the evolution for the calculation with rai/ra2 = 1.5. The capture into the 
9\ 0° and 62 ~ 180° configuration and the passage to the 61 62 ^ 0° configuration are as for 
m\/m2 = 0.3 in Figure 1. However, when the increasing ei reaches ^ 0.17, the libration centers 
depart by tens of degrees from 0°, leading to stable libration of the 9j far from either 0° or 180°. 
Asymmetric librations persist for 0.17 ^ ei < 0.44, while 62 changes from increasing to decreasing 
with increasing ei. The system returns to the configuration with both 61 and 62 librating about 
0° and both ei and 62 increasing smoothly for sufficiently large ei(> 0.44). In Figure 2c we show 
the asymmetric resonance configuration near t/-P2,o = 4 x 10^ in Figure 3. The angle between 
the lines of apsides is slightly more than 90°, and conjunctions occur when neither planet is near 
periapse or apoapse. It should be noted that conjunction means equal true longitudes of the planets, 
which is not the same as equal mean longitudes when 9i and 62 are not 0° or 180°. The evolution 
shown in Figure 3 is characteristic for 0.95 < mi/m2 ^ 2.55, with the minimum (maximum) ei for 
asymmetric librations smaller (larger) for larger mi/m2 (see also Fig. 5a). 
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(a) 01 = 0°, 02= 180° 


(b) 0, = 02 = O° 




{ * / ~r 




mj/mg = 0.3, t/Pgg = 6x10^ 


(c) 


(d) 


^^^^^^^^^^^^^ 


( \ * / 


m/m3= 1.5, t/P2,„ = 4x10' 


m/nij = 5.0, t/Pjo = 2x10^ 



Fig. 2. — Examples of 2:1 resonance configurations that can be reached by differential migration 
of planets with constant masses and initially nearly circular orbits. The dashes on the ellipses 
representing the orbits mark the positions of the periapses, and the planets represented by the 
small dots are shown at conjunction, (a) Anti-symmetric configuration with 9i 0° and 62 ~ 180° 
at small eccentricities. The eccentricities of the orbits are exaggerated so that the positions of the 
periapses are more visible. (6) Symmetric configuration with 9i m 0° near t/P2,o = 6 x 10^ in 
Fig. 1. (c) Asymmetric configuration near t/P2,o = 4x 10^ in Fig. 3. {d) Asymmetric configuration 
with intersecting orbits near i/-P2,o = 2 x 10^ in Fig. 4. 

There appears to be a narrow range in mi/m2 (2.55 < m\/m2 < 2.75) for which the system 
changes from 0\ 0° and O2 ~ 180° to asymmetric librations without passing through the 6\ 
O2 ~ 0° configuration as e\ increases, but the system does return to the 61 82 ^ 0° configuration 
for sufficiently large ei. Figure 4 shows the evolution for the calculation with mi/m2 = 5, which is 
characteristic for m,i/m,2 > 2.75. The system changes from 9i ~ 0° and 62 ~ 180° to asymmetric 
librations and never reaches the 9i k, O2 k, 0° configuration. In the asymmetric libration branch, 
62 increases monotonically with time but e\ does not. The system eventually becomes unstable at 
t/P2,o = 2.06 X 10^. The asymmetric libration configurations in Figure 4 with 62 > 0.35 have the 
additional property that the orbits intersect. An example — the configuration near t/P2,o = 2 x 10^ 
— is shown in Figure 2d. 
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5x10^ 10' 

t/P2.0 



Fig. 3. — Same as Fig. 1, but for mxjm^ = 1.5. The sequence of resonance configurations after 
resonance capture is {61,02) ~ (0°, 180°) (0°,0°) asymmetric librations (0°,0°). 




Fig. 4. — Same as Fig. 1, but for mijvni = 5. The sequence of resonance configurations after 
resonance capture is {61,62) ~ (0°, 180°) — asymmetric hbrations. 
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Fig. 5. — Trajectories of 2:1 resonance configurations from differential migration of planets with 
constant masses and initially nearly circular orbits. The trajectories are labeled by m\/m2- The 
total planetary mass (mi + 1112) /mo = 10~^, and the outer planet is forced to migrate inward 
with 02/02 = — 10~^/P2- (o) The ei-e2 plane in linear scales. The configurations above the dash 
for nil/ 1712 = 3, 5, and 10 and the configurations between the two dashes for mi/ni2 = 2.3 and 
2.65 have intersecting orbits. The diamond, squares and circles on the curves with mi/m2 = 0.3, 
1.5 and 5 indicate the equilibrium eccentricities for eccentricity damping of the outer planet with 
K = 162/621/102/02! = 1, 10 and 100, respectively. See text for K = 1 and mi/m2 = 0.3 and 1.5. 
The configurations indicated by the arrows on the curves with mi/m2 = 2.65 and 3 are used as 
initial conditions for calculations in §3.2 and Figs. 7 and 8. (6) The small eccentricity region of the 
61-62 plane in logarithmic scales. The solid lines show the relationship in eq. (7) from the resonant 
perturbation theory to the lowest order in the eccentricities. The dashed line is the boundary below 
which the 9i = 0° and 62 = 180° configuration exists, as determined by Beauge et al. (2003). 

It should be noted that the symmetries of the problem without forced migration mean that 
every stable 2:1 resonance configuration with centers of libration about 9i and 62 has a counterpart 
with —9i and —02- However, the forced migration causes small offsets in the libration centers, and 
there could be a preference for evolution into one of the asymmetric libration branches, depending 
on the migration rate and the libration amplitudes. As in the examples shown in Figures 3 and 4, 
all of our calculations that yielded asymmetric configurations have 61 becoming negative when 61 
initially departs from 0°. 

A convenient and informative way to summarize the results for different mi/m2 is to plot 
the trajectories in the 61-62 plane (Beauge et al. 2003; Ferraz-Mello et al. 2003). Figure 5 shows 
the trajectories for several values of mi/m2 between 0.1 and 10, with Figure 5ffl in linear scales 
and Figure 56 showing the small eccentricity region in logarithmic scales. The trajectories for 
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mi/m2 < 0.95, 0.95 < mi/m2 < 2.55, 2.55 < mi/m2 < 2.75, and mi/m2 > 2.75 [with the 
sequences of Hbration centers (^i, 62) summarized in Table 1] are clearly distinguished in the 61-62 
plane. The transition from {Oi,02) « (0°,180°) at small eccentricities to (^1,^2) ^ (0°,0°) for 
mi/m2 ^ 2.55 occurs when 61 ~ 0.1 and 62 ~ 0, and the transition to asymmetric librations for 
Tni/m2 > 2.55 occurs before ei reaches 0.1, with the critical value of ei smaller for larger mi/m2 
(see Fig. 56). Beauge et al. (2003) have determined the boundary in 61-62 space below which the 
9i = 0° and 62 = 180° configuration exists, using an analytic model for the resonant Hamiltonian 
(with an expansion in eccentricities). Their result is well approximated by the relationship 62 = 
— 0.36ei + 0.035, which is shown as the dashed line in Figure 56. It is in good agreement with the 
boundary in our numerical results. 

For coplanar orbits, the equation for the variation of the periapse longitude is 



dt mjejy^GmoOj dej 

where the disturbing potential 

$ = -^^^<^(/?, 61,62,^1, ^2), (6) 
0'2 

(f is a function of the indicated variables, and /3 = 01/02 ~ 2~^/^, if we consider only the resonant 
terms and neglect terms of order [(mi + m2)/mo]^ and higher (see, e.g., Lee & Peale 2002; Beauge 
&; Michtchenko 2003). To the lowest order in the eccentricities, dzui/dt = /3ni (777,2/777,0) Ci/61 and 
dw2/dt = — 7i2(777i/mo)C'2/e2 for 61 = 0° and 62 = 180°, where nj are the mean motions and 
Ci{(3) = —1.19 and C2(/?) = +0.43 are, respectively, the coefEcients of the 61 cos ^1 and 62 cos ^2 
terms of the disturbing potential (Lee & Peale 2002). Thus, for 2:1 resonance configurations with 
sufficiently small eccentricities, the precessions of the orbits are retrograde, and the requirement 
that the orbits on average precess at the same rate implies the following relationship between the 
forced eccentricities: 

62/61 = -P'^HC2/Ci){mi/m2) = 0.29(7711/7712). (7) 

Table 1. Sequences of 2:1 Resonance Configurations from Differential Migration of Planets with 

Constant Masses 



Planetary Mass Ratio 




Libration Centers (^1,^2) 


7771/7772 < 0.95 


(0°,180°) 


^ (0°,0°) 


0.95 < 777i/7772 < 2.55 


(0°,180°) 


(0°,0°) asymmetric librations ^ (0°,0°) 


2.55 < 7771/7772 < 2.75 


(0°,180°) 


asymmetric librations — > (0°,0°) 


2.75 < 777l/7772 


(0°,180°) 


— asymmetric librations 
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Our numerical results have finite widths in ei and 62 due to finite libration amplitudes, but they 
are on average in good agreement with equation (7) (solid lines in Fig. 5b) for sufficiently small 
eccentricities. We terminate the solid lines in Figure 56 at ei = 0.12/(mi/m2 + 3), which is an 
estimate of where significant deviation from equation (7) would occur, based on the position of the 
local maximum of 62 for mi/m2 <C 1 and the boundary in 61-62 space for the 9i = 0° and O2 = 180° 
configuration. 

As we just mentioned, the orbital precessions are retrograde for the resonance configurations 
with small eccentricities. The precessions remain retrograde throughout the trajectories shown in 
Figure 5a for mi/m2 > 2.75. On the other hand, all of the trajectories for mi/m2 < 2.75 pass 
through a point with ei = 0.714 and 62 = 0.326, where dwi/dt = dw2/dt = and the precessions 
of the orbits change from retrograde to prograde. (Lee & Peale 2002 have already noted the change 
from retrograde to prograde precessions for systems with masses like those in GJ 876.) According to 
equations (5) and (6), the masses only appear as overall factors in dwj/dt, with dwi/dt oc m2/y^mo 
and dw2/dt oc mi/^mo. Thus for small amplitude librations and nearly constant ej and 9j, the 
combinations of ej and 6j for which dzui/dt = dw2/dt is nonzero depend only on mi/m2, and the 
combinations for which dwi/dt = dw2/dt = arc independent of the planetary and stellar masses. 
Since all of the cases in Figure 5a with mi/m2 < 2.75 have 9i and 62 librating about 0° at large 
61, the only point on the 61-62 plane where the trajectories can intersect is a combination of ei and 
62 for which dtvi/dt = dw2/dt = 0. 

In our discussion above of the case with mi/ 1712 = 5, we pointed out that the asymmetric 
libration configurations with 62 > 0.35 have intersecting orbits (see, e.g.. Fig. 2d). Ferraz-Mello et 
al. (2003) have stated that configurations with the apocentric distance of the inner planet greater 
than the pericentric distance of the outer planet [i.e., ai(l-l-ei) > 02(1 — 62) for mi, 777,2 <S tzq] may 
have intersecting orbits, but this condition is necessary but not sufficient (it is both necessary and 
sufficient only for anti-aligned configurations such as those discussed in §4.2). We have determined 
which of the configurations shown in Figure 5a actually have intersecting orbits. For mi/m2 = 3, 
5, and 10, the configurations above the dash on each trajectory have intersecting orbits. However, 
configurations with intersecting orbits are not limited to mi/m2 > 2.75. For mi/m2 = 2.3 and 
2.65, the configurations between the two dashes on each trajectory also have intersecting orbits. 



3.2. Effects of Migration Rate, Total Planetary Mass, and Eccentricity Damping 

We consider next the effects of a migration rate on the disk viscous timescalc similar to that 
in equation (4) and a larger total planetary mass. We have performed differential migration calcu- 
lations similar to those in §3.1, but with 02/02 = — 10^^/P2 and (mi -|- m2)/mo = 10^'^ and 10^^, 
and the trajectories of stable resonance configurations in the 61-62 plane are shown in Figures 6a 
and 66, respectively. Figures 5a and 6 show that the basic trajectories depend only on mi/m2 and 
not on the total planetary mass, in agreement with the analysis in §3.1 based on the requirement 
of equal precession rates. However, the libration amplitudes and the point at which a sequence 
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Fig. 6. — Same as Fig. 5a, but for 02/02 = — 10~"^/P2 and (a) (mi +m2)/mo = 10~^ and (6) 
(mi +m2)/mo = 10~^. The basic trajectories depend only on mi/m2, but the libration amphtudes 
and the point at which a sequence becomes unstable can be sensitive to the migration rate and the 
total planetary mass. 



becomes unstable can be sensitive to the migration rate and the total planetary mass. 

For (mi +m2)/mo = 10^'^, the faster migration rate (Fig. 6a) leads to larger libration am- 
plitudes when a system with mi/ni2 ^ 0.95 enters asymmetric libration, and the point at which 
a system with mi/m2 > 2.75 becomes unstable is slightly different from that shown in Figure 
5a. All of the trajectories in Figure 66 for (mi + m2)/mo = 10~^ are wider than those in Figure 
5 a because the larger eccentricity variations generated by the larger total planetary mass before 
resonance capture lead to resonance configurations with larger libration amplitudes. On the other 
hand, although the trajectories in Figures 6a and 6b are for the same migration rate, those in 
Figure 66 with mi/m2 > 0.95 do not show any significant increase in the libration amplitudes 
when a system enters asymmetric libration, because the larger total planetary mass means that the 
resonant interaction between the planets is stronger relative to the forced migration. It is not clear 
in Figure 66, but with the use of Jacobi coordinates (see §2), we are able to see the capture into 
the 9i ~ 0° and 62 ~ 180° configuration at small eccentricities when we examine time evolution 
plots (like those in Figs. 1, 3, and 4) for these calculations with (mi + m2)/mo = 10~^. 

In Figure 66, many of the trajectories for mi/m2 < 2.3 show deviations from the trajectories in 
Figure 5a at ei > 0.75 due to rapid increase in the libration amplitudes before the system becomes 
unstable (see §4.1 for a more detailed discussion of a similar case). The cases with mi/m2 > 2.65 
become unstable near the points where the orbits become intersecting (dashes on the trajectories 
in Fig. 5a). Note, however, that the case with mi/m2 = 2.3 passes through the configurations 
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with intersecting orbits without becoming unstable. The fact that a sequence becomes unstable at 

a certain point docs not mean that all configurations beyond that point are necessarily unstable, 
especially for smaller libration amplitudes. We have used the asymmetric configuration, with 
non-intersecting orbits, indicated by the arrow on the curve with mi/m2 = 2.65 in Figure 5a 
as the starting point for a calculation in which mi + m2 is increased from 10~^mo at a rate of 
din (mi + m2) /dt = 10~®/P2,o while m\/m2 is kept constant. It remains stable to the end of the 
calculation when (mi + m2)/mo = 10~^. When we used instead the asymmetric configuration, 
with intersecting orbits, indicated by the arrow on the curve with mi/m2 = 3 in Figure 5o as the 
starting point, the configuration becomes unstable when (mi + 1712) /itiq exceeds 0.0033. 

As we mentioned in §1, either eccentricity damping or the termination of migration due to 
nebula dispersal would lead to a system being left somewhere along a sequence in Figures 5 and 
6. There is significant uncertainty in both the sign and the magnitude of the net effect of planet- 
disk interaction on the orbital eccentricity of the planet because of sensitivity to the distribution 
of disk material near the locations of the Lindblad and corotation resonances (e.g., Goldreich & 
Sari 2003). Nevertheless, hydrodynamic simulations of two planets orbiting inside an outer disk 
have shown eccentricity damping of the outer planet, with K = \e2/e2\/\a2/a2\ ~ 1 (Kley et al. 
2004), while an explanation of the observed eccentricities of the GJ 876 system as equilibrium 
eccentricities requires K ~ 100 (Lee & Peale 2002). To study the effects of eccentricity damping, 
we have repeated the mi/m2 = 0.3, 1.5, and 5 calculations shown in Figures 1, 3, and 4 with 
K = 1, 10, and 100. For K = 10 and 100, the eccentricities reach equilibrium values in all cases, 
and they are indicated by the squares and circles in Figure 5 a. For K = 1, while the eccentricities 
in the mi/m2 = 5 calculation reach equilibrium values indicated by the diamond in Figure 5a, 
those in the mi/m2 = 0.3 and 1.5 calculations are still increasing at the end of the calculations 
when ei « 0.68 and 0.70, respectively. 



4. OTHER RESONANCE CONFIGURATIONS 

4.1. Symmetric and Asymmetric Configurations 

The bifurcation of the trajectories in Figure 5a at mi/m2 ~ 2.75 leaves an empty region in the 
61-62 plane. In this subsection we show that there are stable 2:1 resonance configurations in the 
empty region of Figure 5a that cannot reached by the scenario considered in §3 (i.e., differential 
migration of planets with constant masses and initially nearly circular orbits). 

In Figure 7 we show the results of a calculation in which the asymmetric configuration indicated 
by the arrow on the curve with mi/m2 = 3 in Figure 5a is used as the starting point and mi/m2 
is decreased from 3 to 0.1 at a rate of d\n{mi / 7712) / dt = — 10~^/P2,o while (mi -|- m2)/mo = 10~^ 
is kept constant. As mi/m2 decreases, 61 and 62 remain near the initial values, 62 changes little, 
and ei increases (i.e., the sequence extends to the right in the 61-62 plane). Similarly, we have 
used the asymmetric configuration indicated by the arrow on the curve with mi/m2 = 2.65 in 
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Fig. 7. — Evolution of the eccentricities, ei and 62, mean-motion resonance variables, 9i and 62, 
and mass ratio mi/m2 for a calculation in which the configuration indicated by the arrow on the 
curve with mi/m2 = 3 in Fig. 5a is used as the starting point and mi/ 1712 is decreased at a rate of 
dln{mi / m2) / dt = — 10~^/P2,o- The configurations with mi/m2 = 1, 0.3, and 0.1 indicated by the 
dashed lines are used as initial conditions for calculations in Figs. 9 and 11 and §4.3. 
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Fig. 8. — Same as Fig. 7, but for a calculation in which the configuration indicated by the arrow 
on the curve with mi/m2 = 2.65 in Fig. 5a is used as the starting point and mi/m2 is increased 
at a rate of d\n{rn\ / 1712) / dt = 10~^/P2,o- The configurations with m\/m2 = 3, 5, and 10 indicated 
by the dashed lines are used as initial conditions for calculations in Figs. 10 and 11 and §4.3. 
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Fig. 9. — Resonance configurations from forward and backward migration calculations using the 
mi/m2 = 0.3 configuration from Fig. 7 as initial conditions. The results from the calculations with 
^2/02 = — 10~®/P2 and 10~®/P2 are plotted along the positive and negative time axis, respectively. 
These resonance configurations are different from those in Fig. 1 from differential migration of 
planets with constant masses and initially nearly circular orbits. 



Figure 5a as the starting point for a calculation in which m\/m2 is increased from 2.65 to 10 at a 
rate of d\n{rnx / 1712) / dt = 10~^/P2,o- The results are shown in Figure 8. As mi/m2 increases, 9i 
and 62 remain near the initial values, ei changes little, and 62 increases (i.e., the sequence extends 
upward in the 61-62 plane). Thus we have found for each mi/m2 between 0.1 and 10 a resonance 
configuration distinct from those shown in Figure 5o. 

For a given mi/m2, we can then use the resonance configuration from either Figure 7 (if 
mi/m2 < 2.75) or Figure 8 (if mi/m2 > 2.75) as the starting point for "forward" (a2/«2 < 0) 
and "backward" (02/02 > 0) migration calculations to search for other resonance configurations. 
Figure 9 shows the results from the m\/m2 = 0.3 calculations with 02/02 = — 10'~^/P2 and 10'~^/P2 
along the positive and negative time axis, respectively. In both calculations we have integrated the 
system until it becomes unstable. The asymmetric resonance configurations in Figure 9 are clearly 
different from the symmetric and anti-symmetric configurations in Figure 1. Figure 10 shows the 
results from the mi/m2 = 5 calculations. The calculation with 02/02 = — 10~^/P2 was terminated 
when 61 reaches 0.8, while the calculation with 02/02 = 10~^/P2 was terminated when both 61 and 
62 circulate. We shall not consider further configurations with only 62 librating, but it is interesting 
that those in Figure 10 at —8.5 x 10^ < ^/^2,o ^ —4.6 x 10^ are similar to those in Figure 4 at 
t/P2,o ^ 6 X 10^ with both 61 and 62 librating. The configurations in Figure 10 with both 61 and 62 
librating (at i/P2,o ^ —4.6 x 10^) are either symmetric or asymmetric, and they are clearly different 
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Fig. 10. — Same as Fig. 9, but for calculations using the mi/m-2 = 5 configuration from Fig. 8 
as initial conditions. The configurations at i/-P2,o ^ —4.6 x 10^ with both 9i and 62 librating are 
different from those in Fig. 4. 

from those in Figure 4. 

In Figure 11 we show the trajectories in the 61-62 plane of stable resonance configurations 
(with both 61 and librating) from the forward and backward differential migration calculations 
with initial conditions (triangles in Fig. 11) from Figures 7 and 8 for mi/m2 = 0.1, 1, 3, and 
10. The trajectories for mi/m2 = 0.3 and 5, whose time evolutions are shown in Figures 9 and 
10, respectively, are simply between the appropriate pairs in Figure 11 and are not shown to 
avoid crowding. Unlike the cases with mi/m2 = 3 and 5 (see Fig. 10 for the latter), where 62 
increases monotonically with increasing ei for the Oi ^ 62 ^ 0° configurations at large ei, the 
trajectory for mi/m2 = 10 shows a maximum in 62 at 61 « 0.72. This turns out to be due to 
the libration amplitudes of 61 and 62 becoming too large. As the example in Figure 10 shows, 
the libration amplitudes of 9i and O2 increase as the system is driven deeper into resonance. The 
semi-amplitudes do not exceed 10° in the cases with mi/m2 = 3 and 5, but for mi/m2 = 10, the 
configurations with 61 > 0.68 have semi-amplitudes of Oi and 62 > 10° (although 9i and 62 are 
nearly in phase so that the libration amplitude of 03 = 61 — 62 remains small and the trajectory 
in the ei-62 plane remains narrow). We have repeated the m\/m2 = 10 calculation by starting at 
ei = 0.61 with initial conditions that correspond to smaller libration amplitudes. The resulting 
trajectory is shown as the gray curve in Figure 11, and it does not show the decrease in 62 at large 
ei. 

Like the trajectories in Figure 5a for m\/m2 < 2.75, the trajectories in Figure 11 for mi/m2 = 3 
and 10 (the gray curve for the latter) also pass through the point with ei = 0.714 and 62 = 0.326, 
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Fig. 11. — Trajectories in the ei-e2 plane of 2:1 resonance configurations from forward and back- 
ward migration calculations with initial conditions {triangles, oriented to indicate the direction for 
forward migration) from Figs. 7 and 8 for mi/m2 = 0.1, 1, 3, and 10. See text for the calculation 
that produces the gray curve for m\/m-2, = 10. The resonance configurations are distinct from those 
in Fig. 5a. All of the configurations for mi/m2 = 0.1 and 1 and the configurations to the left of 
the dashes on the trajectories for mi/m2 = 3 and 10 have intersecting orbits. 

where dtui /dt = dw2 /dt = and the precessions of the orbits change from retrograde at smaller 
ei to prograde at larger ei. In Figure 11, all of the configurations for mi/ 1712 = 0.1 and 1 and the 
configurations to the left of the dashes on the trajectories for mi/m2 = 3 and 10 have intersecting 
orbits. 

The configurations shown in Figure 11 clearly occupy a different part of the 61-62 plane com- 
pared to those in Figure 5 a. They are found by a combination of calculations in which mi/m2 
is changed slowly and differential migration calculations, and cannot be reached by differential 
migration of planets with constant masses and initially nearly circular orbits. 

4.2. Anti-symmetric Configurations with 9\ « 180° and 02 ~ 0° 

As we mentioned in §1, Ji et al. (2003) and Hadjidemetriou & Psychoyos (2003) have discovered 
stable 2:1 resonance configurations with 9i ~ 180° and 62 ~ 0° for systems with masses similar 
to those of the HD 82943 system announced initially by the Geneva group [mi/m2 = 0.54, (mi + 
m2)/mo = 2.3 x 10"^]. We have derived from the results of Ji et al. (2003) a 6*1 180° and 02 ~ 0° 
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nii/mj = 0.54, 0^ = 180°, 6^ = 0° 




e, = 0.465, = 0.490 



Fig. 12. — Anti-symmetric resonance configuration with 9i ~ 180°, ^2 ~ 0°) f'lid intersecting orbits 
used as the starting point for calculations in Fig. 13. 
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Fig. 13. — Evolution of the eccentricities, ei and 62, mean-motion resonance variables, 61 and 62-, 
and mass ratio m\/m2 for calculations in which the configuration shown in Fig. 12 is used as the 
starting point and m\/m2 is increased to 10 and decreased to 0.1. The results from the calculations 
with dln[mi/m2)/dt = 10~^/P2,o and — 10~^/P2,o are plotted along the positive and negative time 
axis, respectively. The configurations with mi/m2 = 0.1, 1, and 10 indicated by the dashed lines 
are used as initial conditions for calculations in Fig. 14. 
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Fig. 14. — Trajectories in the 61-62 plane of 2:1 resonance configurations with 61 ~ 180° and 62 ~ 0° 
from forward and backward migration calculations with initial conditions {triangles, oriented to 
indicate the direction for forward migration) from Fig. 13 for mi/m2 = 0.1, 1, and 10. These 
resonance configurations are distinct from those in Figs. 5 a and 11. The trajectories terminate 
above the dashed line, where ai(l + ei) = 02(1 — 62) and the planets would collide at conjunction. 

configuration with mi/m2 = 0.54, (mi + m2)/mo = 10~^, and small libration amplitudes, and it 
is shown in Figure 12. The 9i ~ 180° and 62 ~ 0° configuration has intersecting orbits, with the 
periapses nearly anti-aligned and conjunctions occurring when the inner planet is near apoapse and 
the outer planet is near periapse (and closer to the star than the inner planet). 

We have used the configuration in Figure 12 as the starting point for calculations in which 
mi/m2 is increased to 10 and decreased to 0.1. The results are shown in Figure 13, with the 
calculations with dln{mi/m2)/dt = 10~^/P2,o and — 10~^/P2,o along the positive and negative 
time axis, respectively. The libration centers remain at 9i = 180° and 62 = 0°, and ei decreases 
{62 increases) with increasing mi/m2- We have thus found for each mi/m2 in the range 0.1 < 
mi/ 1712 < 10 a resonance configuration with 9i « 180° and ^2 ~ 0°. As in §4.1, we can then use 
the resonance configuration for a given mi/m2 from Figure 13 as the starting point for forward 
and backward migration calculations to search for other resonance configurations with the same 
m\/m2- Figure 14 shows the trajectories in the 61-62 plane of stable resonance configurations from 
the migration calculations with mi/m2 = 0.1, 1, and 10. The initial conditions (triangles in Fig. 
14) arc indicated by the dashed lines in Figure 13. The calculations with 02/0.2 = — 10^^/P2 were 
terminated when 62 reaches 0.7, while the calculations with 02/02 = 10~^/P2 were terminated when 
the system becomes unstable. All of the configurations shown in Figure 14 have libration centers at 
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01 = 180° and 02 = 0° and retrograde orbital precessions. Along the sequence for a given mi/m2, 
the fractional distance between the planets at conjunction, (ai/a2)(l + ei) — (1 — 62), decreases 
with decreasing 62, and the system becomes unstable at a point above the dashed line in Figure 
14, where ai(l + ei) = 02(1 — 62) and the planets would collide at conjunction. We can see from 
Figure 14 that the configurations with 0i ^ 180° and 02 ~ 0° tend to a sequence with ei = and 
62 above 0.37 in the limit mi/m2 — ^ 00. These are the pericentric libration configurations found by 
Beauge (1994) for the exterior 2:1 resonance in the planar, circular, restricted three-body problem. 

4.3. Effects of Total Planetary Mass 

We consider in this subsection the effects of a larger total planetary mass on the stability of 
the resonance configurations found in §4.1 and §4.2 for (mi + m2)/mQ = 10^^. We have used the 
configurations indicated by the dashed lines in Figure 7 as initial conditions for calculations in which 
mi + 777-2 is increased at a rate of din (7771 +7772) /dt = 10~^/P2,o while 7771/7772 is kept constant. 
These configurations with 7771/7772 = 1, 0.3, and 0.1 become unstable when (7771 + 7772)/777o exceeds 
0.0037, 0.0041, and 0.0043, respectively. Similar calculations with the configurations indicated by 
the dashed lines in Figure 8 with 7771/7772 = 3, 5, and 10 as initial conditions show that these 
configurations are stable for (7771 + 1712) /mo up to at least 10~^. 

For the anti-symmetric configurations with 0i ~ 180° and 02 ~ 0°, the minimum 62 (and 
the minimum fractional distance between the planets at conjunction) at which the configurations 
become unstable increases with increasing total planetary mass. We have used the configurations 
at the top of the trajectories shown in Figure 14 with 62 ~ 0.7 and 7771/7772 = 0.1, 1, and 10 as 
initial conditions for calculations in which mi + m2 is increased, and they remain stable to the end 
of the calculations when (7771 +m2)/mo = 10~^. When we used instead the configurations indicated 
by the triangles in Figure 14 as initial conditions, the cases with mi/m2 = 0.1 and 1 arc stable 
for (7771 + m2)/mo up to at least 10~^, but the case with 7771/7712 = 10 becomes unstable when 
(7771 + m2)/mo exceeds 0.0096. 

It is likely that many (if not all) of the asymmetric configurations shown in Figure 9 for 
7771/7772 = 0.3 (which is close to 7771/7772 = 0.32 of the GJ 876 planets) are unstable for the total mass 
of the GJ 876 planets [(7771 -|-7772)/777o = 0.0095], because the configuration at t = is the 7771/7772 = 
0.3 configuration from Figure 7 examined above and it is unstable for (7771 + m2)/mo > 0.0041. 
Thus for a system with masses like those in GJ 876, there are stable 2:1 resonance configurations 
with (6*1, 6*2) ~ (180°, 0°), in addition to those with (6'i,6l2) ~ (0°,180°) and (0°,0°) that can be 
obtained from differential migration of planets with constant masses. However, the planets are 
sufficiently massive that asymmetric configurations like those in Figure 9 are probably unstable. 
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5. DISCUSSION AND CONCLUSIONS 

We have shown that there is a diversity of 2:1 resonance configurations for planar two- 
planet systems. We began with a series of diflferential migration calculations, with planets having 
constant masses and initially circular orbits, and found the following types of stable 2:1 reso- 
nance configurations: (1) anti-symmetric configurations with the mean-motion resonance variables 
^1 = Ai — 2A2 + tni and ^2 = Ai — 2A2 + W2 librating about 0° and 180°, respectively (as in the 
lo-Europa pair); (2) symmetric configurations with both 61 and 62 librating about 0° (as in the 
GJ 876 system); and (3) asymmetric configurations with Oi and 62 librating about angles far from 
either 0° or 180°. Systems with mi/m2 < 0.95, 0.95 < mi/m2 < 2.55, 2.55 < mi/m2 < 2.75, 
and m\/m2 > 2.75 show different types of evolution (Table 1), and their trajectories in the 61-62 
plane are clearly distinguished (Fig. 5). The basic trajectories depend only on m\/m2-, but the 
libration amplitudes and the point at which a sequence becomes unstable can be sensitive to the 
migration rate and the total planetary mass. Where a system is left along a sequence depends on 
the magnitude of the eccentricity damping or the timing of the termination of migration due to 
nebula dispersal. There are also stable 2:1 resonance configurations with symmetric (^1 ^2 ~ 0°), 
asymmetric, and anti-symmetric (^1 180° and 62 ~ 0°) librations that cannot be reached by 
differential migration of planets with constant masses and initially nearly circular orbits. We have 
analyzed the properties of the resonance configurations found in our study. The asymmetric con- 
figurations with large 62 and the 9i « 180° and 62 ~ 0° configurations have intersecting orbits, 
while the 61 k, O2 0° configurations with e\ > 0.714 have prograde orbital precessions. Figures 
5, 11, and 14 can be used to determine whether an observed system could be near a 2:1 resonance 
configuration with both 9i and 62 librating. 

We now describe several mechanisms by which the resonance configurations summarized in 
Figures 11 and 14 can be reached. Since we found the symmetric and asymmetric configurations 
in §4.1 (Fig. 11) by a combination of calculations in which m\/m2 is changed and differential 
migration calculations, they can be reached by differential migration if the planets continue to 
grow at different rates and mi/m2 changes during migration. In particular, if mi/m2 < 2.75 
when the system is first captured into resonance and mi/m2 increases above 2.75 when ei > 0.17, 
the system should cross over to the configurations with mi/m2 > 2.75 in Figure 11. Similarly, if 
mi/m2 > 2.75 when the system is first captured into resonance and m\/m2 decreases below 2.75 
when 62 > 0.49, the system should cross over to the configurations with mi/m2 < 2.75 in Figure 
11. After the planets open gaps about their orbits and clear the disk material between them, they 
can continue to grow by accreting gas fiowing across their orbits from the inner and outer disks. 
If the inner disk is depleted, the situation with the outer planet growing faster and a decreasing 
mi/ 1712 may be more likely (Kley 2000). 

A mechanism by which the 9i ~ 180° and 62 ^ 0° configuration in §4.2 (Fig. 14) can be 
reached is a migration scenario involving inclination resonances. Thommes & Lissauer (2003) have 
studied the capture into and evolution in 2:1 resonances due to differential migration for orbits 
with small initial inclinations. For the case with 1711/1712 = 1, the system is initially captured 
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into eccentricity- type mean-motion resonances only (with both 0i and 62 Hbrating), and the initial 
evolution after capture is similar to that of the planar system discussed in §3. But when ei 
grows to ~ 0.6, the system also enters into inclination- type mean-motion resonances (with both 
2Ai — 4A2 + 2O1 and 2Ai — 4A2 + 2O2 librating, where Oi^2 are the longitudes of ascending node), 
and the orbital inclinations begin to grow. The system eventually evolves out of the inclination 
resonances, and interestingly, the eccentricity resonances switch from the ^2 ~ 0° configuration 
to the Oi f» 180° and 82 ~ 0° configuration at the same time. Thommes Sz Lissauer (2003) have 
noted that the inclination resonances are not encountered until ei grows to ~ 0.6, which requires 
m,i/m2 to be less than about 2 for planets with constant masses. From the results in §3 (Fig. 
5o), we find that the boundary in mi/m2 for ei to reach ~ 0.6 is more precisely mi/m2 < 2.75. 
However, if we allow mi / m2 to change during migration as described in the previous paragraph, 
it would be possible for e\ to grow to ~ 0.6 even for systems with final m\/m2 > 2.75 (see Fig. 
11). Thus it may be possible to reach the Oi 180° and 62 ~ 0° configuration by evolving through 
inclination resonances for non-coplanar systems with a wide range of final mi / 1112 ■ 

Finally, there is a mechanism for establishing mean-motion resonances in planetary systems 
that does not involve migration due to planet-disk interaction. Levison, Lissauer, & Duncan (1998) 
have found that a significant fraction of systems at the end of a series of simulations starting 
with a large number of planetary embryos has mean-motion resonances. More recently, Adams 
&; Laughlin (2003) have performed a series of iV-body simulations of crowded planetary systems, 
with 10 giant planets initially in the radial range 5 AU < a < 30 AU. After 10^ yr, there are often 
only two or three planets remaining in a system, and a surprisingly large fraction (~ 10%) of the 
systems have a pair of planets near resonance (roughly equally distributed among the 2:1, 3:2, and 
1:1 resonances), although not all of the resonant pairs will survive on longer timescales. If there 
is no additional damping mechanism, the resonance configurations produced by multiple-planet 
scattering in crowded planetary systems are likely to have moderate to large orbital eccentricities 
and relatively large libration amplitudes. A more detailed analysis of the resonances produced by 
this mechanism is needed, but there is no obvious reason to suspect that it cannot produce any of 
the 2:1 resonance configurations with moderate to large orbital eccentricities presented in §3 and 
§4, provided that the configuration is stable for relatively large libration amplitTidcs. 

The number of multiple-planet systems is expected to increase significantly in the next few 
years. Our results show that a wide variety of 2:1 resonance configurations can be expected among 
future discoveries and that their geometry can provide information about the origin of the reso- 
nances. 

It is a pleasure to thank Stan Peale for his contributions to the early stages of this work and for 
many useful discussions. I also thank G. Laughlin and E. W. Thommes for informative discussions. 
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